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RESONANCES AND LOCALIZATION 
IN MULTI-PARTICLE DISORDERED SYSTEMS 

VICTOR CHULAEVSKY 



Abstract. This is a complement to our earlier work [1] where a new eigenvalue 
concentration bound for multi-particle disordered quantum lattice systems was 
obtained. Here we show that the new bound leads to a simplified proof of multi- 
particle spectral and dynamical localization. 



1. Introduction. The model and the motivation for this paper 



We study multi-particle quantum systems in a disordered environment, usually re- 
ferred to as Anderson-type models. Consider an iV-particle tight-binding Hamiltonian 
H Vi uM in the Hilbcrt space £ 2 {Z Nd ), 

N 

H v ,u = ^ + gy + lJ = ^2^A^+gV(x j ,uj)) + U, (1.1) 

i=i 
where A is the nearest-neigbor lattice Laplacian, 

O: AW*(x) = Atf>tt(si,...,a:tf)= Y. *(xi,...,Xj +y,...,x N ), 

CH ! yeZ d :\y\ = l 

V : Z d x — > R is a random field relative to a probability space (f2, J 7 , P), g > is 
a constant measuring the "amplitude" of the potential V, and U is the multiplication 
operator by a function U(x) which we assume bounded, but not necessarily symmetric. 
^^ ' In addition, we assume that U has a finite range ro < oo. In the case of a two-body 

interaction generated by an interaction potential U 1 - 2 ^ : Z d — > R, 

U(xi,...,XAr) = ^2U {2) {\x t ~Xj\), 

X ■ 

this hypothesis is equivalent to the following condition: suppC/' 2 ' C [0, ro]- The as- 
sumptions on the random field V are described in Sect. 11.31 Unless otherwise specified, 
the boldface symbols denote objects relative to multi-particle systems. 

Given any finite cube Ci(u) := {x e Z Nd : ||x — u|| < L}, we will consider a 
finite- volume approximation of the Hamiltonian H 

Hc L (u) = H t^ 2 (Ci(u)) with Dirichlet boundary conditions on 9Cj,(u) 

acting in the finite-dimensional Hilbert space I 2 (Cl(u)). In [5] the following "two- 
volume" version of the Wegner bound was established for pairs of two-particle oper- 
ators Hc t ( u ), He ,( u ') such that L > L' and dist(Cz,(u), Cz/(u')) > 8L: if v is the 
continuity modulus of the marginal distribution function Fy , then 

P{dist(a(H Ci(u) ),a(H Ci , (u0 )<e} < (2X + l) 2d (2L' + l) d v{2e). (W2) 
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2 V. CHULAEVSKY 

1.1. More efficient eigenvalue concentration bounds. In [HIE] the MPMSA was 
used to prove spectral localization (i.e., exponential decay of eigenfunctions) in the 
strong disorder regime. Aizcnman and Warzel [H|3] used the FMM to prove directly 
dynamical localization (hence, spectral localization) in various regions in the parameter 
space, including strong disorder, "extreme" energies and weak interactions. 

However, due to a highly correlated nature of the potential energy in multi-particle 
systems, it was difficult to obtain optimal decay bounds for eigenfunctions in terms of 
some norm in Z Nd . This difficulty had been analyzed by Aizenman and Warzel [2] and 
served as the motivation for a more efficient eigenvalue concentration bound proven in 
[4] (cf. Thm. ITTT1 below) . Here we prove spectral and dynamical localization for multi- 
particle systems with the potential V(x; to) satisfying the hypotheses of Thm. 11.11 for 
\g\ large enough |J (strong disorder regime). 

We would like to emphasize that the novelty of the present paper consists in the 
observation that the proof of dynamical (and spectral) localization for multi-particle 
systems, based upon Thm. 11.11 can be obtained by a minor modification of the single- 
particle version, once the key MSA bounds are established for the multi-particle system 
in question. In the author's opinion, this opens a way to numerous extensions of existing 
"single-particle" techniques to disordered systems with interaction. 

1.2. Basic geometrical definitions. Consider the lattice (Z ) = Z Nd , N > 1. 
Vectors x = (xi, . . . , xn) G % Nd will be identified with TV-particle configurations in Z d . 
We use below the max- norm || • l^ in R nd D Z nd , n > 1: 

||x||oo = ||0ci,---,Zn)lloo := max max \x\'\. 

3e[[i,n]]ie[[i,d]]' 3 

This norm canonically induces the notion of diameter for subsets of M. nd and Z nd , 
denoted below as "diam". We denote by ID) the "principal diagonal" in (Z ) : 

© = (x e Z Nd : x = (x, . . . , x), x € Z d }. 

We will often use the standard notation [[a, b}] := [a, b] (~l Z. 

Taking into account the symmetry of the potential energy V(xi ; w) + • • • + V(xn ; uS) , 
it is natural to introduce also the " symmetrized" distance 

ds(x,y) := min ||x-y||. 

rGtow 

Given a cube Cl(u) = {x : ||x — u|| < L) c Z Nd , we define three kinds of its 
" boundaries" : 

3-C L (u) = {x: ||x-u|| = L}, 

a+C L (u) = {x: ||x-u||=L + l}, 

aC L (u) = {(x,x') : x e d-C L (u),x' e a+C L (u), ||x- x'|| = 1}. 

1.3. The main result on multi-particle eigenvalue concentration. Introduce the 
following notations. Given a parallelepiped Q C Z d , we denote by £q( w ) the sample 
mean of the random field V over Q, 



An adaptation of our method to the case of weak disorder, at "extreme" energies, is the subject 
of a fortgeoming manuscript by T. Ekanga [10] , 
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and introduce the " fluctuations" of V relative to the sample mean, 

rj x = V(x,w) -£q(u), x e Q. 

We denote by $v,Q the sigma-algebra generated by {r) x , x G Q; V(y; ■), y £ Q}, and by 
F^(- \3v,q) the conditional distribution function of £q given $v,Q : 

^(*|3v,0):=PU<3<*l3v,<}}. 

We will assume that the random field V satisfies the following condition: 

(CM(i/)): For any R > there exists a function vr : R+ — > R+ vanishing at and 
such that V Q C Z d with diam(Q) < R the conditional distribution function F^(- |3V,q) 
satisfies 

Vt,sGM, eB8Bap\Fs(t\$v,Q)-Fe(8\!Sv,Q)\<VR(\t-8\). (1.2) 

This condition may prove useless for applications if vr(s) 4- too slowly as s J. 0. 
For this reason, we assume below, in addition to ( 11.21) . that 

VR 



(t) < Const R A hT B \t\- 1 , \t\<l, (1.3) 



for some A < oo and sufficiently large B. Alternatively, a stronger condition of Holder- 
or Lipshitz-continuity can be used: for some A < oo and b > 0, 

v R {t) < Const R A \t\\ \t\ < 1. (1.4) 

Note that in the particular case of a Gaussian IID field V with zero mean and unit 
variance, £q is a Gaussian random variable with variance \Q\ , independent of the 
"fluctuations" n x , so that its probability density exists and is bounded: 

p (Q (s) = \Q\^(2n)-^e-^ < \Q\^(2n)~^, 
although ||pf Q ||oo grows with \Q\, and so does the continuity modulus of F^ Q . 

Theorem 1.1 (Cf. @]). Let V : Z d x n ^ R be a random field satisfying (CM(i/)). 
Then for any pair of N -particle operators He ,(u')> He »(u")> < L',L" < L, satis- 
fying dg(u', u") > 2NL, and any s > £/ie following bound holds: 

P{dist(a(H Ci , K) ),a(H Ci „ („»))) < s} < |C L '(u')| • |Cl»(u")I"l(2*). (1-5) 

An adaptation of our method to a class of bounded potentials V(-;u>), for which 
(CM(i/)) does not hold, is the subject of a forthcoming manuscript by M. Gaume |12) . 

2. MULTI-PARTICLE MSA 

2.1. Decay of resolvents in finite cubes: the key bound. 

Definition 2.1. Given a sample H(w) of an iV-particle Hamiltonian of the form (|1.1[) . 
a cube Cl(u) is said (E, m)-non-singular ((E, m)-NS), with £el and m > 0, if 

max max |(H Cl(u) - ^(x, y)l < e ^(™^), ( 2 .1) 

x: ||x— iiWKL 1 /" y£o~C L (\i) 

where 

7 (m, L, n) := mt(l + Zr 1 /*)"-"* 1 , l<n<N. 
Otherwise, it is called (E, m)- singular ((E,m)-S). 
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Traditionally, the fc-dependence is put in the decay exponent, m = rrik, which is 
re-calculated recursively. We prefer to keep fixed the decay parameter to, while the 
actual decay bound depends on Lfc and on N explicitly, through the function 7. 

The main result of the multi-particle MSA, used for the derivation of spectral and 
dynamical localization, can be formulated as follows: 

(DS.I,k, n): For any pair of 2nLk-distant n-particle cubes Cl,.(x), C_L fc (y) the 
following bound holds: 

¥{3EeI: C Lk (x) and C Lk (y) arc (E, ro)-S }<V 2 "" +1 . (2.2) 

2.2. Geometric resolvent inequalities. The second resolvent identity implies the 
following formula for the resolvent Gc i ( u )(x,y; C) 0I the lattice Schrodinger operator 
in a cube Cl(u) D Cf(w), L > I + 1: for any complex £ which is not an eigenvalue of 
operators H Ci(u) and H c<(u) , and x e Q(w), y G C L (u) \ Q(w), 

G Ci( u)(x,y;C)= E G Ci(w) (x,v;C)G Ci (u)(v',y;C) (2.3) 

(v,v')e3C,(w) 

yielding the well-known Geometric Resolvent Inequality (GRI): 



|G Ci (u)(x,y;C)| <\dCi{w)\ max |G c , (w) (x, v; C)| |G Cl(u) (v',y;C)| 

(v,v')E()C((w) 

< \dC e (w)\ ( max |G c ,( w )(x, v;C)| j , max |G Cl (u)(v',y; C)|- 

\v£<)-C((w) J v'£3 + C((w) 



(2.4) 



Furthermore, if \& n is an eigenfunction of the operator H Cl / u \, then for any complex 
£ which is not an eigenvalue of operator H Cf ( u ) the following identity holds true (we 
will call it the Geometric Resolvent Inequality for eigenfunctions) : 

*„(x)= J2 G C£(w) (x,v;C)*„(v'), x€C/(w), (2.5) 

(vy)eac,(w) 

yielding the GRI for eigenfunctions: 



|¥ n (x)| < |0C,(w)| max |G c , (w) (x, v; C)| max |*„(v')|. (2.6) 

\vea-c ( (w) / v'ea+c f (w) 

The derivation of the GRI can be found, e.g., in the review (T4"] . 

2.3. Base of the induction in N. In the course of the induction in n, we will assume 
that (|2.2p is established for all n < N — 1 and all k > 0, and then prove it for n = N, 
at each scale L&. The parameter p > will be required to be sufficiently large (cf. 
(O), and the exponent 2 N ~ n+1 p in the RHS of (|2T2"j) shows that the strongest bound 
is required for the single-particle systems (n = 1). It follows from the results of single- 
particle MSA that p can be made arbitrarily large, if the disorder parameter g is large 
enough. The single-particle FMM provides even an exponential bound of the LHS in 
(|2.2[) for n — 1. See, e.g., the articles [T1I51HTJ. the review [2] and the monograph |15) . 
As n grows, the upper bound in (|2.2[) becomes weaker, hence, simpler to prove. 



2.4. Base of the scale induction, for a fixed N. 

Lemma 2.1. There exists g* < 00 and positive functions m*{g),p*{g) defined for 
Iff I > 9* such that 

(1) m*(g),p*(g) ->■ +00 as \g\ ->■ 00; 
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(2) the property (DS 1,0, N) holds for the operators H (JV) (w) = A (JV) + .oV(cj) + U 
wit/i 771 = ?7i* (<j>) and p = p* (g) . 

This statement can be proven for A-particlc Hamiltonians in the same way as for 
N = 1; see, e.g., [8llll|. In fact, the non-random interaction U (and/or a non-random 
component of V) can simply be ignored in the proof. Note for a reader not familiar with 
the conventional MSA techniques that m*(g) ~ 0(ln \g\), for \g\ 3> 1. The asymptotic 
behavior of p*(g) depends upon the regularity of the marginal distribution function Fy 
of the random field V. 

2.5. "Radial descent" bound for resolvents. In this subsection we state an an- 
alytic result which does not rely upon single- or multi-particle structure or random 
nature of the potential energy of the Hamiltonian H. All constants involved depend 
only upon the dimension of the lattice. In order to emphasize this, we change here our 
notations and do not use boldface symbols. Since V is arbitrary here, in applications 
to the multi-particle systems we can assume that V contains also an inter-particle in- 
teraction energy. In a certain sense, Lcmma l2. 21 below encapsulates an argument going 
back to [8l[TT] and used since then in many papers. 

Definition 2.2. A bounded function / : A -> C on a subset A C Z", n > 1, is called 
{£, q, S, c)-subharmonic, with £ £ N, q > 0, «S C A, c > 1, if for any x € A \ S such that 
Ct(x) C A, one has 



while for x G S 



|/(x)| <q max |/(y)|, 

\\ y -x\\=e 



I f(x)\ < a max I f(v)\ 



Lemma 2.2. Consider an (£ 7 q,S,c)-subharmonic function f onCi,(u) C Z n . Suppose 
that the c£-neighborhood of the set S can be covered by a collection A of annuli 

Ai = C bi (u)\C ai (u), 

of total width W A . Then for any r £ [W A +£,L-W A +£] 

max |/0)| < q [(L-r-W(A))/£]-l max |j( j|_ 
xeCV(u) y£C L (u) 

In particular, 

|/(0)| < q \(L-w(A))/£\-i max | /(y ^ 
yec L (u) 

The proof, based on a "reverse induction" in r ("radial descent") is given in Sect. O 
A direct application of this bound to the resolvents, making use of the GRI, leads to 
the following statement. 

Lemma 2.3. Consider a lattice Schrodinger operator A + V (x) in a cube Cj, k+1 {u) C 
Z n , n > 1. Fix an energy E E I and suppose that the cL^ -neighborhood of all (E,m)- 
singular cubes of radius L& inside C'L k+1 {u) can be covered by a collection A of annuli 
Ai = Cbi(u) \ C ai (u) of total width W A . Suppose also that CL k+1 {u) does not contain 
any E -resonant cube of radius L > Lk (including itself). Then, for any fixed value of 
the constant c and Lq > Lq(c, n) large enough, the cube C_L fc+1 (u) is (E,m)-NS. 

We will see below (cf. Sect. I2.7J) that, in application to MPMSA, with sufficiently 
high probability a cube of radius Lk+i does not contain more than 4 singular cubes 
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of radius Lfe which are pairwise 2iVX/ £ -distant, in which case all singular cubes can be 
covered by a collection of annuli of total width O(Lk). 

We stress that Lemma 12.31 is purely " deterministic" and does nor rely upon singlc- 
or multi-particle structure of the potential. 

2.6. Localization bounds for decomposable systems. We need a simple property 
of quantum systems decomposed in a union of perfectly non-interacting subsystems. 
Loosely speaking, it says that if two distant subsystems u' G Z™ , u" G Z n are 
"localized" and do not interact, then their union u = (u',u") is also "localized". 

Definition 2.3. 

(1) Let n' € {1, . . . , N- 1}, k > and u' = (u x , . . . , u n >) G Z"' d . Given a bounded 
interval Icl and m > 0, the n'-particle cube C}™ (u') is said (m, 7)-tunneling 
((m, I)-T, for short) if 3 E G I and there are 2./VLfc_i-distant n-particle cubes 

C^vxj.cg^to) C c£V) which are (E,m)-S. 

(2) An .ZV-particle cube C^ (u) is said (to, Impartially tunneling ((m, J)-PT) if, 
for some permutation t G &n (acting on the components of the vectors u = 
(ui, . . . , un)) and some n', n" > 1, it admits a representation 

C W(r(u))=c£ i (u')xC^ i (u") 

u' = (ui,...,u n >), u" = («n'+i,.-.,«iv), where cither C^^^u') orC^^u") 
is (m, 7)-T. Otherwise, it is said (m, 7)-NPT. 

Lemma 2.4 (Cf. Lemma 3 in [7]). Fix an interval / C I and an energy E G J. 

Consider an N -particle cube C L (u) = C)™ (u') x C)™ (u"), wii/i n',n" > 1, and a 
sample of the potential V(-;u)) such that 

(a) p (UC^ (u'),UC { ^p (u")) > 2L k + r ; (r = the range of interaction) 

(b) C^ (u) is E -non-resonant; 

(c) C^u') and Ci r 2(u") are(m,I)-NT; 

(d) i/ie inductive assumptions {(DS.I,k', n), A:' > 0, n < iV — 1} hold true. 

Then the cube C L (u) is (E,m) -non- singular. 



Lemma 2.5 (Cf. Lemma 5 in [7]). If Lq is large enough, then 

1 r _ p2 N-(«-l) + l +2A , da 

2 



{C Lk (u) is (m,I)-PT} < ± L - P *— + ™^ (2.7) 



Indeed, the "tunneling" property says that one of the projection cubes C[ (u') 
contains at least two distant singular cubes of radius L/t_i, so that Lemma 12.51 follows 
from the inductive assumptions {(DS.I, k',n), k' > 0, 1 < n < N — 1} combined with 
a simple bound on the number of pairs of points in C)™ (u') by i|C^(u)| 2 . 

Introduce the following random variables, depending upon the samples V(-;uS) : 

K pl (CL k+1 (u), J) = the maximal number of 2iVLfc-distant PI cubes of radius L k 
in Ch k+1 (u) which are simultaneously [E, m)S, with E E I, 

K F 7 (Ci fc+1 (u),7) = the maximal number of 2iVLfc-distant FI cubes of radius Lk 
in Cz,j. +1 (u) which are simultaneously [E, m)-S, with E E I. 



'{Jf"(C %1 (u),i)>4} < ^r +iNd < U- k H (2.9) 
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The maximal collections of resonant and distant FI (or PI) cubes are, of course, not 
uniquely defined, but their (maximal) cardinality is well-defined, since there is only a 
finite number of choices for cubes Cl,.(v'^), . . . , ClJv^) inside Ci t+1 (u). 

Corollary 1. If p > Nda/{2 - a) and l^- a )- Nda > 2 , then 

P { K PI (C Lk+1 (u), I) > 2 } < ^ L -<*+™** < 1 L'H (2.8) 

2.7. Pairs of non-decomposable (FI) cubes. 

Lemma 2.6 (Cf. [7]). I/p > 2Nda/(2 - a), then 

1 
4[-/c+r 

For the proof, it suffices to notice that distant FI cubes give rise to independent 
samples of the potential V. 

As was mentioned before, we see that with sufficiently high probability K PI +K FI < 
3 + 1 = 4, and Lemma l2.3l can be used to prove non-singularity of a non- resonant cube 
of radius Lk+i with at most K FI + K PI < 4 pairwise 2iVXfc-distant singular cubes of 
radius L^ inside it. Indeed, 0(Lfc)-neighborhood of all such cubes can be covered by a 
collection of annuli of total width O(Lk). 

Note that, with a = 3/2, it suffices to require that p > 6Nd and Lq > 2; then the 
hypotheses on p in Corollary [T] and in Lemma 12.61 arc automatically satisfied. 

Corollary 2. Assume that the bound (DS.I,k, N)(Eqn (|2.2|) holds for some fc > 0. 
Then the bound of the form (DS.I, k + 1, N) also holds for all 2N L^ + i~ distant pairs 
of FI cubes. 

Proof. Consider two 27VXfc + i-distant FI cubes Ci fc+1 (x), CL fc+1 (y). Both of them arc 
(E, m)-S for some E € I only if at least one of the following events occurs: 

(1) for some Eel, both Ci t+1 (x) and Ci fc+1 (y) are E-R; 

(2) K FI (C Lk+1 (x),I) > 4 or K PI (C Ljb+1 (xj,/) > 2; 

(3) K FI (C Lk+1 (y),I) > 4 or K PI (C Lk+1 (y),I) > 2. 

By Thin. 11.11 the probability of the event (1) can be bounded by \L~^.J' 1 (actually, by 
any power of L^ +l , provided that (|1.3[) is fulfilled with B > large enough). In the 

case where (|1.4|) holds, one obtains even a stronger bound by e~ k + 1 , (3' > 0. Next, it 
follows from Corollary [5] and Corollary Q] that the probability of the event (2) (as well 
of the event (3)) is bounded by 

I t ~ 2 p -l — t - 2 p ^ It - 2 p 

Therefore, the cubes Cl h1 (x) and C_L fc+1 (y) are simultaneously (E, m)-singular, for 
some E S /, with probability < ( j + |) £fe + i < ^fc+ii as required. □ 

2.8. Pairs of decomposable (PI) cubes and mixed pairs. 



Lemma 2.7. Assume that the bound (DS.I,k, N)(Eqn (|2.2[) holds for some k > 0. 
Then the bound of the form (DS.I, k + 1, N) also holds for all 2N Lk+i-distant pairs 
of cubes Ci t+1 (x) and Ci fc+1 (y) ; of which at least one is "partially interactive". 
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Proof. Without loss of generality, assume that C^ fc+1 (y) is PI. The cubes Ci fc+1 (x) 
and Ci fc+1 (y) are simultaneously (E, m)-singular, for some E € I, only if at least one 
of the following events occurs: 

(1) for some E € I, both Ci H1 (x) and Ci fe+1 (y) are E-R; 

(2) C Lt+1 (x) is FI, and K FI (C Lk+1 (x),I) > 4 or K PI (C Lk+1 ( X ),I) > 2; 

(3) the PI cube Ci fc+1 (y) is partially tunneling. 

Naturally, the option (2) is to be considered only for a mixed pair of cubes: if both of 
them are PI, it suffices to analyze only the cube CL fc+1 (y). 

The probabilities of the events (1) and (2) can be estimated exactly as in the proof 
of Corollary [21 so that their sum does not exceed \L~ kJ ^ l . By inductive assumptions 
{(DS.I,k', n), k' > 0} on systems with n < A — 1 particles, the event (3) has proba- 
bility bounded by jL^ + ^. Combining these estimates, the lemma follows. □ 

Now the key bound of the MPMSA, (DS.I, k, N), is established for all k > 0. 

3. From the MSA bounds to multi-particle localization 

It has been known since more than twenty years that the principal MSA bounds - 
in the single-particle theory - imply the spectral localization; see the original papers 
[8|lllj. In fact, as was pointed out in 6,7, the main argument used in such a derivation 
is not specific to single- or multi-particle structure of the random potential. Note for a 
reader familiar with the traditional, single-particle version of Lemma 13 . 1 1 given below, 
that the main argument used in the single-particle context applies to multi-particle 
systems, with one minor technical modification: instead of disjoint pairs of cubes, one 
should consider 2 A L^- distant pairs of cubes at each scale L^. 

In essence, owing to the new eigenvalue concentration bound given by Thm. 11.11 the 
derivation of 2-particle spectral localization from MSA bounds of the form (DS.I, k, N) 
described in [6] extends - in a fairly simple way - to any A > 2. 

Lemma 3.1 (Cf. [7]). Fix an interval Icl and suppose that the bound (DS.I, k, N) 
(cf. Eqn (|2.2[) ) holds for all k > 0. Then with probability one, the Hamiltonian Yi\ N > (ui) 
has pure point spectrum, in I, and all its eigenf unctions ^ n {uj) with eigenvalues E n G / 
decay exponentially fast at infinity: 

3C„(w) : Vx€Z JVd |*»(x;w)| < C„(w) e~ m||x|1 . 

Proof. As it is well-known, spectrally-a.e. generalized eigenfunction >& of a lattice 
Schrodinger operator H is polynomially bounded, 

|*(x)| < Const ||x|| A , A < oo, 

so it suffices to show that every polynomially bounded solution of the equation H\l/ = 
E^ with E G / is exponentially decaying at infinity. 

If \l/ ^ 0, then there is a point u G Z Nd where \P(u) ^ 0. We start the analysis of 
the decay properties of ^ by finding the smallest integer kg > such that the cube 
Gh k (0) contains all points r(u), r G &n, obtained by permutations of the components 
Uj of the vector u = (ui, . . . , ujv)- 

The first observation is that for some k\ > ko and all k > k\ the cubes Ci k (0) 
must (E,m)-S. Indeed, assume otherwise. Then there exist an infinite number of 
values of k (hence, arbitrarily large values of Lfe) such that Cl & (0) is (JS,m)-NS, 
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with u G Ci fc _ 1 (0). Then the non-singularity property implies that for all points 
x G Ci, s ._ 1 (0), including x = u, one has 

|*(x)|<e^( m ' L - JV |c>C Lfc (0)| max |*(y)| 

yea+Ci, fc (o) 

<e- mL »0(L[ N - 1)d+A ) —> 0. 



Lfc-S-OG 

This would mean that ^(u) = 0, which is impossible by our choice of the point u. 
Next, set a' = 9/8 < a and for all k>k\ consider the events 

£k = \ 3 A G / : C 2LQ ' (0) contains two 2iVXfc-distant (A, m)-S cubes of radius Lk > 

and 

Q:= |J f|(°\^)- 

k>k± j~>k 

Since ¥{£ k } = 0(L 2 k Ndaa '~ 2p ) andp > 3Nda > Ndaa', it follows from Borcl-Cantclli 

lemma that P< $7 > = 1, and for every u> G fl (i.e., with probability one) there exists 

k% > k\ such that no cube Cz, +1 (0) with j > k2 contains a pair of (E, m)-S cubes of 
radius Lj at distance < 2NLj. 

Further, introduce a sequence of annuli 

A fc = C ir+i (0)\C L a,(0), fc>0, 
and let x G Z, Nd \ C L <*' (0) = Ufc>fc 2 Afc. (Note that it is not necessary to estimate 

k 2 

|M/(x)| for x G C L „/ (0) otherwise than by a constant.) Suppose that x G Aj for some 

k 2 

j. Observe that the cubes Ci fc (0) and Cl (0) arc 2iVLj -distant and u G 17, so that one 
of them is (£,m)-NS. Since fc > fc 2 , C Lj (0) is (£,m)-S, hence, C Lj (x) is (£,m)-NS. 
Moreover, all cubes of radius Lk outside CL k +2NL k (0) are also (E, m)-NS. This implies 
that the function *f? is (Lk , q)-subharmonic in the cube Cr(x), with 

„ _ e -m(L k +0(lnL k )) < e ~ m ( L & + 2 L ' /4 ) ^ 

if ifc is large enough, and 

i? := ||x|| - (L k + 2NL k )), ||x|| > L 9 k /S » L k . 

It is easy to see that i?/||x|| = 1 — 0(L k ), and applying Lemma \2. 31 we obtain, for 
L k large enough, 

as required. D 

Recall that the validity of the hypothesis of Thm. B~T1 is established in Sect. 12.31 We 
come, therefore, to our main result on the A^-particle spectral localization: 

Theorem 3.1. Consider the random operators H> '(w) of the form (|1.1|) . Suppose 
that U is bounded and V satisfies CM(v). Then there exists g* < oo such that if 
\s\ > 9* ; then with probability one the operator H(w) has pure point spectrum, and all 
its eigenf unctions *& n (u>) with eigenvalues E n G / decay exponentially fast at infinity: 

3C7„M :Vx6Z M |*n(x;w)| < C n (uj)e' m ^. 
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Remark 3.1. While the statement of Thm. 13. H is similar to that of the main result of 
[7], a detailed analysis shows that the actual bound on the random constants C n (w) 
obtained in [7] depends upon the position of the localization center x n for the corre- 
sponding cigcnfunction \&„(u;); Thm. [TTT1 rules out such a dependence. 

4. Strong dynamical localization 

To prove the dynamical localization, in addition to the assumption (|1.3|) on the 
random potential V{-\ u) we need the following hypothesis: for any finite interval Jcl 

3k = K (I,N,d) <oo: P{tr(Pj(H Cx ,(u)(w))) >CL« Nd } < L~ B ' , (4.1) 

where -Fj(-) is the spectral projection on / and B' > will be required below to be 
sufficiently large, depending on other parameters. It is readily seen that the trace of 
P/(rIc [ ( u )) grows not faster than linearly in |Cl(u)|, if the random potential V is 
bounded from below. Since we would like to allow also Gaussian potentials for which 
CM(i/) becomes most simple, we allow a faster rate of growth. Observe that 

tr (H Cl ,(u)M) = tr (A Cl(u) + U Ci(u) ) + g tr (V Cx ,( u )M), 

and 

TV 

tr(Vo L („)(w))= J2 J2 V ^'^- 

xec t (u)i=i 
Therefore, for a large class of marginal distributions including Gaussian ones, the re- 
quired bound follows from standard results for the sums of IID random variables. 

4.1. Strong dynamical localization. Here we follow closely the scheme which is 
well-described in [T5]; see also [H, 13 . Specifically, Propositions 1-6 below correspond 
to assertions at Steps 1- 6 from Sect. 3.4 in [15] . 

Introduce the operator X of multiplication by max- norm in £ 2 (Z Nd ): 

(X*)(x) := ||x|| *(x), x e Z Nd . 

Theorem 4.1. Assume that the property (DS.I, k, N) holds true for a given N > 1 
and for all k > 0, with p > (3Nda + as)/2, s > 0. Then the random operators 
H( N > (w) feature strong dynamical localization in the energy interval I: for any finite 
subset K C 1 Nd and any bounded measurable function r\ with supp n C /, 

E[||X s 7 ? ( J tfH)l K ||]<oo. 
Note that the condition p > 2Nda/(2 — a) was required to prove (DS.I, k, N). 
Since the validity of the hypothesis of Thm. l4~T1 is established in Sect. 12.31 it implies 
the iV-particle dynamical localization for the Hamiltonians H(u>): 

Theorem 4.2. Assume that the random field V : 7L d x O — > R satisfies the condition 
CM(i/) and |j (|4.ip with B 1 > 2p — Nda. Assume also that the initial scale estimate 
(DS.I, 0, N) is fulfilled for some interval 7 C K with p > (3Nda + as)/2, s > 0. 
Then the random operators U^ N ' (ui) feature strong dynamical localization in the energy 
interval I: for any finite subset K C Z Nd and any j] G L°°(R) with supp n C I , 

E[\\X S n(H(u)) Ik || ]< co. 

In particular, H> ' (w) features complete dynamical localization with any given value 
of s > 0, if g is large enough: \g\ > g*(s), g*(s) < oo. 



Recall that I I4.1I I is not required for random potentials bounded from below. 
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Now we will describe the strategy of the proof of Thm. 14.11 Our main goal here 
is to show that the proof, making use of the simpler and more general eigenvalue 
concentration bound (|1.5|) . is very close to that used in the single-particle context. 

4.1.1. "Bad" and "good" events. Fix s > 0. We will always assume that 

p > max {2Nda/(2 - a), (3Nda + as)/2} (4.2) 

and set b = b(p, N, d, a) := 2p — 2Nda. For each j > 1 consider the events 

Sj = {u: 3E £ I 3y,z £ C 4{N+1)Lj+1 (0) such that d s {y,z) > 2NLj 
and C Lj (y),C Lj . (z) are (m, £)-singular}. 

and (in the case where the random potential V is not bounded from below) 

Tj = {w : tr (Pi(H 0ii+a(u) )) > CL$i }. 

Further, for k > 1 denote 

and consider the annuli 

M fc = C 4(A r +1 ) Lfc + 1 (0) \ C 4 (Ar + i) Lfc (0). 

Proposition 1. Under the assumption (|4.ip iw£/i B' > 2p — 2Nda 
Vfc > 1 P{^ bad) } < c(a,d,p,N)L~ {2p ~ 2Nda) . 

Proof. The number of pairs y,z£ C4(jv+i)l +i (0) figuring in the definition of the event 
Sj is bounded by |(4(JV + l)L j+1 ) 2Nd = C(N,d)L 2 j Nda , and (DS.IJ.N) says that 

P { C Lj (y), C Lj (z) are (m, £)-singular } < LJ 2p , 

while for the event Tj we have the bound (|4.1[) . Now we require that the exponent S' 
be large enough, so that P { Tj } < LJ p+ a and 

p | n<£ ad) } < E ( p { s j ■ } + p { % }) ^ E ConstL j Nda ■ 2 L i 



-2p 



r~ 6 



i>fe 


j>k 


E^ 6a "' 


< Const L-} 2p ~ 2Nda) 


j>fe 





(4.3) 



D 

4.1.2. Centers of localization. By Thm. l3~Tl there exists a subset $7i c f2 with P { tti } — 
1 such that for any u £ f2i the spectrum of H^^w) in / is pure point. Fix welli 
and let <&„(w) be a normalized eigenfunction of H(w), with eigenvalue E n (uj) £ I. We 
call a center of localization for $„ every point x n (w) £ Z Nd such that 

|*„(x„(w))| = max |*„(y)|. (4.4) 

y£Z Nd 

Since 3? n G ^ 2 (Z 7Vd ), such centers always exist and, due to the normalization ||$„|| = 1, 
the number of centers of localization x„ ia for a given n must be finite. 

Proposition 2. There exists fco such that for all w S fii anrf k > kg, if one of the 
centers of localization x„ ja /or an eigenfunction \& n ™£/i eigenvalue E n £ I belongs to 
a cube Ci t (x), t/ien t/ie cw&e Cj, t+1 (x) is (m,E n )-S. 



12 V. CHULAEVSKY 



Proof. Assume otherwise. Then the GRI for the eigenfunctions (|2.6p combined with 
the non-singularity condition (|2.1[) implies that, for L& large enough, 



|*n(x„, a )| < e-^ m ' L »+^ |0C Lfc+1 (x)| max |*„(y)| 

< max |* n (y)|, 

y eS+c i/o+1 (x) 

which contradicts the definition of a center of localization. □ 

We will denote by x Mj i is the center of localization closest to the origin; again, it 

might be not unique, but this is irrelevant for the proofs. Fix ko as in Prop. [2] and set 

n { k 900d) = n 1 \n<£ ad) . 

Proposition 3. There exists jo = jo(m, a, d) large enough such that for j > jo, j > k 

andx„ tl G C Lj+1 (0) 

K 1 ~ 1 c < „ +1)1 . i+a (!o))*n\\<i- 

Proof. Using the annuli M&, we can write 

ii(i-ic 4(N+1) , j+1 (o))* n |i 2 = E II im,*^ E E i*»(y)i 2 

i>j+2 j>j+2 ygM, 

Fix i > j + 2. The cube Cl^^O) 3 Cz, +1 (0) contains the center of localization 
x Mj i, so, by Prop. [3J Ci 4 (0) is (i?„,m)-singular. By construction of the event Q}? 00 , 
k < j, the cube C^y) must be (i?„,m)-NS. In turn, this implies by the GRI for 
eigenfunctions that |<I?„(y)| 2 < e ~ 2mLi . Now the claim follows from a polynomial (in 
L t ) bound on the number of terms in the sum X^vpm (')• ^ 

Proposition 4. There exists C4 = 04(771, d, k) such that for uj G ^tf°° and j > k the 
following bound holds: 

#{n: x„£C lj+1 (0)}<c 4 I-f. (4.5) 

Proof. Observe that we have 

E (lc^(o)ft(H)lo lHa (o)*», *») <tr (lc Lj+2 ( 0)J P/(H)). (46) 

x„.iec Lj+1 (o) 

It suffices to show that each term in the LHS is bigger than i. Indeed, by Prop. OH 

(lc^ +a (0)ft(H)l Oli+a (0)*»,*n) 

= fl C ,. +2 (0)P/(H)*„, * n ) - (l Ci . +2 (0)P/(H)(l-lc t . +2 (0))*n, *r 

> lc L , + ,(0) *„, *n) - 3 = (*», *n) - ((1 - 1 Cl , + 2 (0)) *n, #„) - 4 



4 
> I. 



□ 



4.1.3. Eigenfunction correlator bounds. 



Proposition 5. There exists an integer k\ = k\(K,Lo) such that for all k > k\, 

7 7 (^(a;))(x,0)|< e ^^-/ 2 ||7 ? || oo . (4.7) 



u G n { k 9 ° od) and x G M fc 



RESONANCES AND LOCALIZATION IN MULTI-PARTICLE SYSTEMS 



13 



Proof. Without loss of generality, assume that \\t]\\oo ¥" 0: 

|| ?? || 00 |- 1 ?? ( J ffH)(x,o)|< J2 INUr 1 !^)!!*^*)!!*^)! 

n: E„el 
< Yl |*n(x)ll*n(0)| + X; E l*n(x)||*„(0)|. 



(4.8) 



n:E n £l 
x„,!6M fe (0) 



j>k n:E n eI 

x„ i i£Mj(0) 



For k > ko and Lk large enough, the first sum in the RHS can be bounded as follows: 
J2 |#„(x)| |* n (0)| < const L^ d e~ mL - < X -e~ mL -l\ 



n: E n el 
x„,ieM fe (0) 



(4.9) 



since one of the cubes Ct t (x), Ci fc (0) must be (E n ,m)-NS: indeed, these cubes are 
2iVLfc-distant. Next, fix any j > k + 1 and consider the sum with x„ i 6 Mj(0). The 
cubes 0^.(0) and C^(x n] i) are 27VXj-distant and by Prop. [2J for k (hence, j) large 
enough the cube Cj^x^i) is (E n ,m)-S, so for w £ f2?£\° the cube 0^(0) must be 
(E n ,m)-NS. Therefore, 

|*„(0)| <Conste- mL i. 
Using again Prop. [4j we get, for Lk large enough, 



J2 E l*n(x)||*„(0)| < Const J2e- mL] L« K Nd < h 



-mL fc /2 



j>k n:E n eI 

x„, 1 gm j (o) 



j>k 



Now the claim follows from (|4.9|) and (|4.10p . 

Proposition 6. Fix k\ as in Prop. [21 Then for any k > k\ and x G Mfc, 



(4-10) 



U 



E 



|| lc £fcW »/(#(")) lo ifc (o) II ] <I^Hoo(cL- 2p+2JVdQ + e- m ^/ 2 ). (4.11) 



Proof. Using Prop. [5] and Prop. [T] we can write 



E 



< 



lc Lfc (x)??(#H) lc lfe (0) II 

lj-^ad) ||l Cifc ( x )77(-ffM) lc Lfc (o)||J +E^l^ SOO d) ||l 0x>k(x )7y(H(w)) lc ifc (o)|| 



E 



a 



4.1.4. Conclusion. Fix a set A" C Z^ and find fc > fci such that A" C ClJO). Then 
E [ ||X S V (H(u)) 1 K || ] < c Nd L s k + ]T E [||X S 1 M , r?(# (")) 1* II] 

<c(fc) + E c ^+i E E[|lc Ii( w)#(w))l 

j>fe \wEMj 

i + ]T if £f da ( L ; 2p+2Nda + e - mL >' 2 ) 



L c tfc (o) 



< CO, 



since 2p — 3Nda — as > 0, and Lj = (Lo) a grow fast enough. 
This completes the proof of Thm. 14.11 
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5. Appendix. Proof of the Radial Descent bound (Lemma | 

Proof. It suffices to consider non- negative functions; otherwise, we replace /by |/|. 

Let a function / : Cx(0) — > M. + be (£, q, S, c)-subharmonic. Introduce "spheres" 
S r = {x : ||x|| = r} and the sets 

cl 

S' = {x : S M n 5 ^ 0}, 5" = |J |J S r+j 

and also 

K = {r >0: C r (0) cC L (0)\5"}. 
Note that if ||x|| = r G 7?., then 

f(x)<q max f(y)<q max /(j/). 

Moreover, for any u with ||u|| G [r — ci", r] the definition of the sets 1Z and 5" implies 
that u G" S, so that 

f(u)<q max f(y)<q max /(y) 

!/: ||l/-u||<* V:||l/||<IMI+< eg ]\ 

- ? ,™ x /(y)- 

j/: ||j/j|<r+£ 

Further, for any u with ||u|| < r — cl there exists a value |f| i? G {^, (1 + c)£} such that 

(5.2) 



f(u) < <Z max /(«) < o max /(y) 

y. ||tf-u||<Jl y. ||l/||<(r-e/)+(l+c)£ 

< q max f{y)- 

y: \\y\\<r+l 



Combining (|5.1[) with (|5.2[) . we conclude that for any r €1Z 

max f(u) < q max f(v). (5.3) 

Construct a sequence of points {r n > 0, < n < M} by recursion: 

ro = L; r n = maxjr G 71 : r < r n -i — £}, 1 < n < M, 

with some M = M(1Z). Set, formally, rjvf+i = 0. We see that, as long as r n > 0, either 
fn = T n -\ —£, or J n := [r n , r n _i — £] C 7Z C . The total length of all non-empty intervals 
J„ is bounded by the total width W(S") of annuli covering S" . Now we can write 

M 

L = y^(7"n_i - r„) 

n=0 
= (fM - »"Af+l) + X! ( r "- 1 ~~ r ") + ^#i n : r « = r «-l ~ ^ 

< £ + W{S") +£#{n : r„ = r n _ a - £} 



yielding 



»r „r „-, L-W(S") 1 i-W^ 

M > #{n : r n = r n _ x - £} > -^- - 1 > -^ - 1. 



Further, W(<S") < Wa, since the annuli Ai G ,4, by assumption, cover the c£- 
neighborhood of the set S. Therefore, 

/(0)< 9 M ||/||oo. 



3 _R = (1 + c)£ iiueS, and R = £, otherwise. 
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More generally, if we stop the construction of the sequence {r n } at n = M' as soon as 
rM'+i < t, for some given r > Wj, + £ > W(S") + I, then we obtain 

f(r)< q M '\\ f \u M >> L - r ; WA -i. 



This completes the proof of Lemma 12.21 D 
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